  %S为S矩阵，J为病种数，I为时间段数，p为p矩阵。
%Qi,m表示i时间段来m个病人的概率，S是对病人的日程安排，这些病人来或不来的概率遵循p矩阵，且相互独立。
%因此i时间段来的病人的最大值根据S的第i行的和确定，我们把范围放大到接到的call数n，这会使Q矩阵比其所需放的数据量大，但无伤大雅。
function Q1=findQ(S,J,I,p,n)
        Q1=zeros(I,n+1);%多给一列来保存m=0的情况。
        for m=0:n%避免重复计算PIE矩阵，先m后i。
        PIE=findPIE(J,m);
        PIE0=PIE;%储存所得结果
          for i=1:I
          PIE=realPIE(PIE0,S,i);%realPIE是一个根据Sij删PIE某些行的程序。
          Q1(i,m+1)=function_Q(i,J,S,p,PIE);
          end
        end
end
